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THE LIFT DISTRIBUTION OF SWEPT -BACK WINGS* 
By J. We is singer 


SYNOPSIS ■ 


Two procedures for calculating the lift distribution 
along the span are given, in which better account is 
taken of the distribution of circulation over the area 
than in the Prandtl lifting-line theory. The methods are 
also applicable to wings with sweepback. Calculated 
results according to two methods agree excellently. Using 
the second more simple method, one needs about 3 hours- 
for the calculation of the lift distribution of a straight ■ 
wing, and about 8 hours for this calculation for the 
swept-back wing. The results are compared with those of 
the Multhopp ‘method and with experimental results. Finally, 
there Is a* brief note on the swept-back wing In sideslip. . 
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INTRODUCTION 


In the present report two methods for determining 
the lift distribution along the span are described, by 
means of which the influence of the wing plan form is 
taken into account better than by previous lifting-line 
theory. Since the new methods are also applicable to 
wings with sweepback, they represent an extension of 
previous theory, at least so far as scope of applicability 
is concerned. 

The method described in the first section is based 
on the assumption of a lifting surface, and therefore will 
be designated as the ’’lifting- surface method,” or 
”F-method.” For a straight rectangular wing, the amount 
of calcxilation required is not materially greater than 
for known methods already in use; it gives a noticeably 
dc Q 

smaller xL than the older methods, and this decrease 

da 

of lift curve slope increases with decreasing aspect 
ratio A and amounts to approximately 8 percent for 
A = 5 • Practically, this fact will be especially important 
for unsymmetrical lift distributions; for instance, one 
can deduce from it a noteworthy decrease of the rolling 
moment due to sideslip resulting from dihedral. In the 
general case of trapezoidal wings with and without sweep- 
back, the required amount of calculation is quite consider- 
able, and consequently one would only use the method In 
special instances for the control of the results from the 
more simple approximate method (L-method). 

The method of the second section is based on a 
slightly modified model of previous lifting line theory, 
and hence will be designated as the ’’lifting-line method” 
or the ”L-m&thod,” In spite of the radical simplification 
of the basic geometrical model compared to that of the 
F-method ( and, consequently, in spite of the marked reduction 
of the required amount of calculation for trapezoidal and 
swept-back wings), the results show an excellent agreement 
with those of the. first method. 

In the third section the results obtained with the 
new methods (the calculations have been subjected to 
numerous checks, but have not been narried'out twice 
independently) are discussed, and compared with experiment. 
In addition, a comparison is made with the Multhopp method 



MCA TM No. 1120 


5 


of* calculating swept-back wings (10), The note on the 
polling moment due to sideslip of swept-back wings 
contained in the fourth section is actually outside the 
scope of the present report, which, except for this not-e, 
is concerned with symmetrical flow incidence, but is made 
here in order to quickly remove a widespread misconception 

In this interim report, intended to make available 
to practice as quickly as possible the results obtained; 1 
to date, the F-metbod is only described in detail for the 
straight rectangular wing, and the detailed application of 
It for general wings is not giv3n, especially 3ince these 
mathematical details are of less interest to the practical 
aerodynamicist , These matters together with some 
supplementary material will be included in a 'later report. 


NOTATION 


coordinates in vortex plane 
(See fig. 1(a).) 


dimensionless coordinates 


2 

A . _ b_ aspect ratio 
F 

b 

\(y) = — =r- local aspect ratio 

t(y) 

Z taper ratio (ratio of root chord to tip chord) 

q> sweepback angle (measured at .i chord line) 

k 

a angle of attack 


x, y or g , 
x ~ 1 72 ’ ^ “ F % 


7 ~ bfe * ^ 


F 

b 

t 


wing area 
wing span 
wing chord 
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£ angle of sideslip 

v incident flow velocity 

c a lift coefficient (a = o a | v2f) 

c m pitching-moment coefficient i M = c m ft v F 

c T rollin 

JU 

a distance of the center of lift of a wing half from 
the plane of symmetry of the wing referred to 
the semispan. 

Y(x,y) circulation density of the hound vortex 
r(y) circulation distribution along the span 

G(y) = dimensionless circulation distribution 

b v 

F(l), F(2.) influence functions for F-method 
L(l) influence function for L~mBthod 

l ~ l(y - rj) argument of influence functions 

I. THE LIFTING SURFACE METHOD (F-METEiOD) 


g- moment 


coefficient 


(l = 


v 2 F 




As remarked in the introduction, this method will 
only be derived here for the simplest case of the straight 
rectangular wing. If 3uch a wing is replaced by a plane 
system of vortices (see fig. 1(a)), in the sense of 
customary lifting surface theory (see for example 
Blenk (2)), and if one denotes the density of circulation 
of the bound vortices by y=(£, n) the induced velocity 
of the point x,y may be calculated as 


W A (x,y) 



P t/2 

^b/2 

i, -£) 2 + (y - n) 2 

1 

J 

-t/2 J 

-b/2 

y - n (x-g ) (y - n) 


6n 


<3£dn /r 

y 


( 1 ) 
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By taking Into account the two equations • . • 

r ^ ^ . 

Y(|_,ri)d£, =r(' r )) — circulation at a 'wing section f’Z) 


r A 

J-t/2 

pb/2 


I vty - 'n )^ dv(c,r,) 

J -b/2 T - TJ 

one can write instead 


dr) = 2 Y (g,y) 


( 3 ) 


— : 

ph/2 ' 

p t/2 

1-1 IS 

II 

•k 

£ 

£ 

r ’ (^)) , i 

tdLal d£ 

X - £ 1=5 

0 

-b/2 ^ 

-t/2 


r fc /2 p\ 


CtT 

VJ 


Lfc/a vi 


»t/2 


■t/2 


r i) 2 . )£Zz T )) 2 . fta(Lni dg d^ 

(4) -• 


(* - §) (y - n) 


6^ 


Y (^,n) is to "be determined by the requirement that 
W^(x,y) ’shall be equal to the prescribed normal comppnent 
V n (x,y) of the incident flow velocity at the wing. If one- 

only takes into account the part enclosed in frame V< A (x,y) 
the condition is- ■ ' .! 


1" 

Sr 



0-t/2 


21 LLal 
x -£ 


V n (x,y) 



pb/2 

hJLSjoJl 


*b/2 


y -rj 


<in (5) 


This is the equation of the two-dimensional problem for' the 
profile section located' - . at the station! . y -having the normal- 
component distribution given on the right hand side of the 
equation. If V n (x,y) is independent of x one may 

obtain from the two-dimensional theory the solution 
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y(£,>7) 



r (y) = trt 


V n (y) 



lb/2 

-b/2 


r f (n) 

7 - rj 


dr) 


( 6 ) 

( 7 ) 


in which the second relation is obviously identical with 
the Prandtl lifting line equation for the determination 
of F(y). According to Pistoles! (11) this relation is 
also valid approximately for normal -component distributions 
which are dependent on x. If for V,~ (y ) one substitutes 
the va?uue of the normal component at the three-quarter 
point (Multhopp (10) has pointed out the importance of 
this idea;, which seems to have been, almost universally 
forgotten, for numerous applications), the approximation 
is rigorously correct if the dependence on x is linear, 
that is to 3ay, the same as it is for a circular-arc 
profile. 


Since the term of Wa(x, 7) not enclosed in the frame 
vanishes for A -4 <» (for, as is proved in (lJ|) the 
equation for the lifting surface goes into equation (5) 
for A-> co ) t and also since it is known that the Prandtl 
equation is in good agreement with experiment down 
to A = 1+, this part may be regarded as a correction term 
to the Prandtl equation, and since it is a correction 
term it need be taken into account only approximately. 

In order to arrive' at an equation for the determination 
of r(n) which will be similar to the Prandtl equation, 
one must first of all get rid of the dependence on x 
in the term in question, and this is done by Pistolesi*s 
approximation by replacing It by Its value at the three- 

quarter chord point (that is, x = ^-), and in the second 

place y(^ti) must be prescribed as a function of g 
for which the most suitable approximation is equation (6). 
One may easily convince himself that the same result will 
be obtained by incorporating both of these procedures in 
the initial equation (1). The following equation is then 
obtained: 
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t//2 v{x -s) g + (y - -n) 2 ' 


-b/2 J-b/2 


(x -g) (y - rj) 




r'(T]>dgdTi (8) 


in whi.ch. x is to be set equal to ^ . 

By introducing dimensionless quantities. 


= S- v. n » id r, . x = £ x 


7 ~ 7 > "H ~ *ir "H* 


§ = 


r(-n) '= b v o(tj), a .= 


we obtain 


a = £ r T ^I dtj 

t-TT -it — 


J-/ " n 


/(x - £ ) + A 2 (y - r/ ! /l -g, 


(x - g)(y -..rj V 1 +£ 


G f (rj - ) d£ dr) 


where. x.' = 0.5 • If for abbreviation, we set 

'v# Q /? - =rfl = i P 1 Vf* - I ) 2 + & (y - 5) 2 /C 


x - £ 


1 + l 
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F 0 * = F~|A(y - Tjff =| (arc si n 5c + ) (12) 

L ~‘ri+y 

p\K{J - rol = - 3 Z-EsH - F(l), l = A(y - (13) 

- ‘-J A(y ~ n) 


equation 10 becomes 


a 


1 + Fp* 
2tt 


n i 




»i 


ff. j - t ii; dt) + 
y - n ' 


p i 


A 

2Tr 


vl-i 



(t)) drj (ll+) 


where F 0 * = 0.88if.7 corresponding to x = 0.5* 

Since F(l) is continuous everywhere, the numerical 
evaluation of the last integral involves no 'fund ardent al 
difficulty. The function F(l) may be calculated with 
the help of elliptical integrals, but practically one would 
get the answer more quickly using numerical approximation 
methods. F(Z-) is an antisymmetric al function whose 
positive branch is shown in figure 2. 

The solution of the integral differential equation ( 1 I 4 .) 
will be carried out in analogy to the Multhopp procedure (9k 
familiarity with which is here assumed, especially familiarity 
with the Multhopp notation (the Multhopp dimensionless 
circulation v is designated here by G-), To this end, we 
need a mechanical integration formula, the proof of which 
is analogous to that of the formula used by Multhopp (9), (10). 

It is : 

If f (rj) is a polynomial ti of degree (2 m + 1), 
that is, of the form 

2H+1 2M+1 

f (r) = > C ri U = y c cos (15) 

v =0 v v -6"" v 

then we have without any approximation: 
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or writing ii/ = cos ^ n 


mr 


f(t)d-,V = 3-^ 


do 

where 


f (^o ) • + f (Vv+n \ JL 

— - g- • + g f(* n ) 


, * 

\{/_ = n — 

^ M + 1 


According to the Multhopp substitution 
2 .m. m 

G(^) = — ----- X- G n fa sin sin uty 


(17) 


(18) 


and using 


2 


f n (t) = sin ^ n co3,(ii|;. 


(19) 


one obtains for G T (ii/) the form. 


G- 1 (^ ) — O n ^n^^ 


( 20 ) 


If this expression Is substituted in. the second integral 
of (lJ|) one gets the following approximation -from 'the 
mechanical integration formula: 


-\ 1 

1 

Zrr 

u-i 



n)_Jo r Cn ) 



m _ 

IlZ 8 un 
n=l 


G. 


n 


( 21 ) 


where g Un = 


2W 


“1 IT 

. p[a(^ - rf)J f n ty) a y 

1—^. ( f™ V+ f »,»*! P v, »*1 + “fc V (ja) 
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in which for the sake of simplicity we have written 


■^np, ^j,) * V 


V 


''At 


arc cos (2p) 

r* 


By taking into account the Multhopp quadrature formula 
for the first integral, one now obtains from equation (lk) 
the following system of linear equations for the 
determination of the 


EL. . 

»/' «v = «u + s b m* <V 0 = X ' 2 ’ • 


m 


( 21 +) 


where 


bp'"" = (1 + Po ,4 )^uv + ^un*" ^ ^ + F o ,4 ^un ~ ^^un 


For symmetrical circulation distributions these equations 
are further simplified to 
m+1 


V a v = %* Sri 


B un G n’ » = x ’ 2 ’ • • • 


m+1 


( 26 ) 


where 


B * 
u 

and 

g 


= (1 + Fo*)t> u „ +AS ot , B, m * = (1 + Fo i! )B,„ - AS 


-1 


uu 

M-l 

~2“ 


DT) 


un 


un 


Un 2 (M + 1) t^ = ° 


F (F 

x np. v up 


p- _) , ‘ u, 
up/ * ^ 


= M + 1 - p. 


(27) 

(28) 


The *np, which appear here are formed from the f n[X and 
are given in table I for m = 7 and M =7# 15* 51 and in 
table II for tn ^ 15 and M = 15*' Their calculation was 
not, made with formulas (19)* but, in particular for the 
case m = M in a considerably simpler manner to which, 
however, no further reference will be made here. 

The method of calculation of the lift distribution of 
a straight rectangular wing with symmetrical angle of attack 
distribution is, then, for the case m ~ 7 as follows: 
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’According. --to the degree of accuracy required, one . 
chooses a certain M, multiplies the differences 
(cos 9 V - cos f ound in table III with the. aspect 

ratio A. then' reads for the points. 7.^^ = A (cos By, - cos 0^,) 
the values P-yu, the function -F(-2>) from figure 2 1 
then forms the difference P-yjj, - and with them using 

H+l 

table I the product sums f n (Fy^ - Py^), which can ; 

jjt =o 

be done on a calculating machine without very much trouble. 
By multiplication with__ -A/2(m + 1) one obtains from 
them the quantities A g un and from them the Multhopp 
coefficients B un , b uv according to formula (27) and also 

the coefficients.' By* - , By n * the system of linear 

equations (26) which may be solved by the Multhopp iteration 
procedure. In doing this it is to be noted that half of 
the coefficients do not vanish as in the case of the 
Multhopp calculations . ■ 

Ixi the case M = m for which most of the l ■ and 
consequently also the P^ are equal but of opposite sign, 

one needs for m — 7 about 3 horn’s for the calculation 
of a lift distribution. For constant, angle of attack over 
the span, the accuracy obtained by setting^ M = m is 
always sufficient, at any rate for .0 % A = 10 although 
the- quantities S un do not corns out very accurately. 

They are, however, partly too large' and partly too small, 

■ so that the lift distribution is hardly affected: for 
instance, for A = 5* a = .1, m = 7 in the calculation 
to 3 decimals there was no distinction between the cases . 

M = 7, M = 15* and M = 31* The increase of accuracy 
resulting frcm the choice of a larger m correspond^ .to 
that of the. Multhopp mothod. 

•The Pistolesi. approximation is rigorously correct for 
normal- component distribution of the form V n ~ o 0 + c^ cos cp 

'where cp = cos"^ x. If the third term of the Pourieir 
development is . also to be considered, the value to be 
taken for ‘Yn is not the value at the three-quarter chord 
point, but the mean of the values at the center of the 
profile and at the trailing edge.- In order to introduce 
this approximation into the above calculation, one must 
regard x as not yet nailed down to a definite numerical 

1 In ordex -1 to improve the accuracy cf reading the values, 
a curve with doubled scale was used. ' 
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T^e course of the calculation_is the sapje as^for equation 
( llj- ) if one replaces F by F • and ' Fo" by S’ o'* = 0.8l8j 
everywhere. The function F(J) is likewise plotted in 
figure 2 . 

A comparison made for a case with A = J shows that 
taking into account the second Pistoles! approximation 
as well as the first does not give anything, at least for 
the degree of accuracy employed in the calculation! c a 

changes by O.J percent. Since, however, after W(l) is 
once calculated, the amount of calculation in the two cases 
is the same; the - ■ systematic calculations for the straight 
rectangular wing were carried out with the function F 
( F-method ). For a number of aspect ratios between 0 and 
10 the lift distributions were calculated on the one hand 
by the Prandtl lifting -line equation, and on the other 
hand by equation (JO); three examples are given in 
figure J. For large A a difference is perceptible only 
at the wing tip where the influence of induction is 
greatest; with decreasing A the difference becomes 
greater, and is in evidence over more and more of the span. 
The na\ir distribution always lies under the old, since 
induction now comes into play more strongly. In addition, 
the corresponding lift coefficients were determined and 
compared with the old values. As A— *0 the ratio 

- — - — tends toward the value O.BS. But oven in the 
1 + F 0 
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range of normal aspect ratios, the deviations are always 
noticeable, as may be seen from figure For purposes 
of comparison, some points for the elliptical wing 
calculated by Helmbold are also included (the distinction 
between rectangle and ellipse for very small A which 
is revealed here is caused by the difference of the basic 
vortex systems) . According to the L-method discussed in II 
we also have for the rectangular wing the limit 0.5 as 
A — ?0« From this figure one can also take a correction 
for the old formulas used to recalculate to a different 
aspect ratio, which correction to be sure ordinarily does 
not amount to very much. 

The new method will give significant differences for 
antisymmetrical angle of attack distributions. Although, 
calculations for such cases are not available, certain 
conclusions may be drawn from figure Ij., For an anti- 
symmetrical lift distribution the lift of a wing-half will 
behave with respect to the induction effect approximately 
like a wing' with half the span; that is, the lift of a 
wing-half calculated as formerly and the corresponding 
rolling, moment are to be multiplied by the ratio 

n t 

O oTp A 

• - ■ read from figure. 5 at the point A • Thus the 
c’ap 

rolling moment due to sideslip caused by dihedral for a 
rectangular wing with aspect ratio ~ 5 would have to be 
about 15 percent less than previously calculated, and 
actually Moeller (8) measured a moment 18 percent too 
small according to the then existing theory. 

VJhen the method, is extended to swept-back wings, a 

number of new difficulties arises, which brings about a 

complication of the formulas and with It an increase of 

the. computational work. For this reason, the explanation 

of the general procedure will not be given in this 

inter in report . For this calculation the basis would be 

the vortex system of the lifting surface as It is shown 

in the example of figure 1(b). The difficulty caused by 

the induced velocity becoming Infinite at the center* of 

the wing is overcome by splitting 'up ; that is, tbe 

circulation distribution G(r>) is substituted in the 
« r ; v 

form Gr(r\) = ET(rf) + c| with a suitable 

constant C or, expressing this in physical terms, the 
singularity of the induced velocity caused by the bound 
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vortices is neutralized by the free vortices. In this 
way we also accomplish something_else at the same time, 
namely that the distribution Q(r\) which is what we are 
solving for, and which is approximated by a Multhopp 
Fourier expression, looks like the distribution of the 
straight wing, so that it can be determined quite 
accurately, even for large sweep angles, with m = ”• 

(See fig. 5 *) The calculation was carried out for 
rectangular wings, A = 5 with sweep angles cp = 0°, 15 0 , I4.5 ; 
the result is shown in figure 6. The purpose of these 
curves is to serve as a basis of comparison with the 
lifting line method which will now be described. 


II. THE LIFTING-LINE METHOD (L- METHOD) 


Given a wing with a straight one-quarter chord line, 
think of the usual model of the lifting line so situated 
at the wing that the lifting line is situated at the one- 
quarter chord line (see fig. 7( a )) and then determine 
(with reference to Plstolesi’s approximation) the 
circulation distribution, so that the vertical component 
of the induced velocity due to this vortex system at the 
thre6-quarter chord line is equal and opposite to the 
corresponding component of the incident flow (to my 
knowledge, this model was first used by Wieghardt (16), 
page 261/262 In a special form) . One has a right to 
expect -that the influence of the free vortices Is pretty 
well. taken into account by this simplified model, since, 
although they are shorter than the ones on the lifting 
surface, they are of constant density and do not decrease 
to zero, as on the lifting surface. Aith reference to 
the bound vortices, it will now be helpful that in the 
case of the infinitely wide plane plate the downwash of 
the lifting surface and the downwash of the lifting line 
located at the one-quarter chord position are equal to each 
other at the three-quarter chord distance . (This result 
is taken from a work by Helmbold which Is not yet 
published.) 


At the general point (x,y) in the plane of the 
vortex sheet the lifting-line model induces the following 
down draft velocity: 


W A (x,y ) = 


1 

5fr 


b/2 


I-b/2 


y - n 


I+37 Vx 2 + (y - t)) 2 


r Kr\) dt) (31) 
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and If one Introduces here the dimensionless notation 


of (9) and also writes A 
becomes 


b 

t(y) 


and 


x 


_ t(y) 
' 2 


this 



u 


-1 



1 + 


+ x^y - - n) 2 


G’(T]) drj 


( 52 ) 


Introducing the notation 

L [>(F - rfO = ^ +*■*<? - ! 

X (y - ri) 


( 55 ) 


we obtain the following equation for the determination 
of G(r})s 


a = 


2 

2fr 


m 


/-I 


i-ial dF + 

y - in 


x 

2tr 


ni 


My - n) 


G f (r,) dT) (3U 


1-1 ~ 


This equation, obviously has the same form as equation (lij.) 
or (30) for the liftings surface method. Consequently the 
numerical treatment Is according to the same scheme as 
given on pages lO/ll, except, that there P Q -”- is replaced 

by 1, F is replaced by L and A by Xy= . The 

antisymmetric function L(Z-) is to be taken from figure 8. 

Ihe application of the model to swept-back wings 
makes no difficulty (compare fig> 7(b)), since there 
.are no ..singularities of the induced velocity at the three- 
quarter chord line. For the down-draft velocity at a 
point (x,y) we have for y = 0 
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w A (x,y) - 


U 


^ b/2 , 

r (n) 

x 4. X - Ini tan cp 

I 

,0 

1 

_ y/{x -Jn tan cp) 2 + (y - n 


d-n 


r> o 


+ & 


3?<*n) ~7 


x + y tan cp 


J-b/2 -]/(x - [ r)| tan cp) 2 + (y - n) 2 


drj 


[*/2 


w 


r(Tj) 


x - y tan cp 


0 


V'x i'll tan cp) 2 + (y - t)) 2 ^ 


dT) (35) 


The last two integrals represent the induced velocity of 
the lifting line, and after integrating them by parts one 
obtains 


r 


Wa = 


n 


l + 


\/(x - jril J :an -C j /; 2 _f (y - r\Y 
x + y tah cp 




2 tan co yen) an 

x^ - y^ tan^ cp! 



'b/2 

— 

1 

] 4- - Ini tan cp ) 2 + (y - n) 2 

tj?r 

ty 

0 

x - y tan cp _ 


IlLnl dn 

y - ti 


(36) 


If now one again chooses the point at which the downwash 
is computed at the three-quarter chord distance, or in 

other words, x = y tan cp + and if one introduces 

dimensionless notation, and also the function 
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one obtains the following equation for the determination of 
G(ti) 




(y.'i) o'(n) an 


<J 8) 


which may be_solved by the same method as the earl-ier equations 
since I>cp(y,r)) is continuous. Compared with the case of 
no sweepfeaok, the^computational work is more tedious because 
the. function Lcp(y,ri) is not a universal function of 
X(y_- T}) as in the case cp = 0°; instead the values 
I‘<p(yu» T )iJ,) must be calculated afresh for every A , t,cp. If 
one writes L^fy,^) as a function of X(y - rj) one gets 

a different function for each point at which th6 downwash 
is computed* y* for Instance if m = J , four different 
functions. In figure 6 the results of calculations by the 
F-method and by the L-method are given for a rectangular 
wing of aspect ratio A= 5 with sweep angles cp = 0 s , 15°, 

'Hie. agreement is very good. If, on the basis of these 
examples, which to be sure should perhaps be increased, one 
assumes a general agreement of the two methods, one may in 
the future use only the less laborious L-method. In contract 
to the Multhopp method (compare section III and also 
figures 11, 12) the convergence of the L-method is also very 
good for large sweep angles; an Increase in the number of 
points from m = 7 to m - 15 does not effect any essential 
change, as may be seen from figure 6. For the same reason 
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ss In the case of the straight rectangular wing (see p. 11) 
one may also calculate gu'in ’-for the pointed and swept -back 
wing by taking- ~M “ m points. Using m = M = 7 one needs 
about 8 hours for the calculation- of the lift distribution 
of the swept-back wing. • • 


III. NUMERICAL RESULTS COMPARISON WITH THE MUDTHOPP METHOD 
AND WITH THE RESULTS .OF EXPERIMENT ' 


In figures 9 and 10 the results of the calculation for 
a rectangular wing of A = 5- and for a trapezoidal wing 
of A = 5» and Z = 2 and different sweep angles are plotted 
For comparison the curves calculated by the Hulthopp method 
(10) ‘for cp = 0° and cp = L\5° and with K = 1 and m = 7 
are Included. The large difference in total lift is 
particularly apparent. This may be explained as follows: - 


Multhopp assumed that the factor of proportionality 
between circulation and angle of attack was independent of 
angle of 3weepback, because experiment showed, at least for 
sweep angles which were not too great, no effect of sweep 
on total lift. The present calculations, ’.however, were 
primarily intended to give as good an approximation as 
possible to the rigorous theory of the lifting surface, in 
order that a solid foundation might be, obtained for the 
estimation of various secondary effects such as boundary 
layer, tip vortices, etc. theoretically, however, sweepback 
must cause a decrease of lift, as one may easily Convince 
himself, and Indeed by a .’factor of cos cp for the wing of 
infinite aspect ratio, while for finite aspect ratio as a 
result of the vortex sheet the decrease is not quite so large 
Based on an approximate' calculation which will not be given 
In detail here, the factor turns out to be ’ 


- A + 2 


‘ %■ 


A 


COS 


<P 


c a’ 


°a’ 


+ 2 
dc t 


a 
da 

for cp 
calculated. 


1 - j!~- x tA-A — in table 5 .the values 
4 - A + j 2 


together with their percentage’ deviations from 

= 0 are given for the examples 'which, were 
The deviations are given very-well by 'the 


r 1 » 

expression * - 'lik . ■ 

2 * A + 2 
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As a second difference between the two methods, one 
sees from' figures 9 and 10 that according to the Multhopp 
distribution the lift is displaced a little more strongly 
toward the wing tips. Table VI gives the position of the 
center of lift on the span according to the Multhopp 

formula ( 1|_. 7 ), (10). Prom the same table one may take the 
difference Aa of the centers of lift for the sweptlback 
wing and the straight wing, the magnitude of which is 
decisive for the question of whether the neutral point of 
the swept-bacic wing may or may not be calculated from the 
lift distribution of the straight wing. ■ Referred to 

average 'wing chord £ the error is tancp Aa ^ 

1 t 2 1 

a quantity which is likewise given in table VI. If a 
maximum error of 1 percent of the average chord is 
permissible, then at least theoretically one must use the 
lift distribution of the swept-back wing beyond cp ~ 20° • One 
sees, moreover, that for cp = If. 5° the difference between 
the Multhopp method, and the lifting-line method is 
considerably greater than the permissible amount. For 
large sweep angles the error introduced by the integration 
formula (ij..7), (10) also plays a role. In order to form 
an estimate of its magnitude, certain values of a were 
determined by planimetry of the corresponding- integral 
areas, and are likewise given in table VI. Here again 
one gets ■ deviations which are too , large for cp = If. 5°. 

Tn an attempt tq explain the differences between the 
Multhopp method and the L-method, the convergence behavior 
of the Multhopp method was examined in certain numerical 
cases. Figure 11 shows one. example (rectangle A = 5> C P = kd°)» 
Multhopp himself points out that, on account of the 
divergence of his integral for .vi^ 'there would be no point 
in the case' of K = .l to increase the number of points 
.m to more than 7* Elis Is confirmed by calculation. Ihe 
distribution calculated for m = 15 is very different from 
the one calculated for m = 7» Also ihen the correction 
function . K (fig. 1 (10)) is used, there Is a marked 
difference at -the center of the wing when m = 7 and 
m = 15* (In this .connection it may be said that all the 
equations of the F-method and of the L-method lead to' the 
form • 
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for suitable HCy,^) from which, using the Multhopp method 
of integration, the following system of linear equations 
is obtained). 


m 

a u = ^uu H uu G u " . 

n=l 



b vn 3uii G n 


• U = 1,2, . ., m 


This method is formally very simple, but converges fright- 
fully slows m = 51 Is not even sufficient. If for the 
swept-back wing one calculates the term due to the sweep 
angle by this method, one has a method which is very 
similar to the Multhopp method for swept-back wing-^with 
K =■ 1, and also has the same convergence behavior. 1 Very 
peculiarly, over the rest of the wing the difference 
between the most exact calculation (m = 15 , K = 1) and 
the roughest calculation (m = 7, K = l)is not excessively 
great. A similar result may be seen for the trapezoidal 
wing with A = 5, Z = 2-, and cp = 1^5°. (See fig. 12.) 
According to this, there does not seem to be much point in 
including the correction factor K In the • Multhopp 
me tho d . 


The question of agreement between theory and 
experiment Is difficult to answer, because at the present 
time there are not many measurements of swept-back wings 
available, and the accuracy of the ones which are avail- 
able Is not always sufficient. The following Is based on 
the rectangular wing measurements by Blenk (5) 

(A = 5 5 cp - 0°, 15 °, 50° ) and Hansen (ij.) (A = I|., 8 ; 
cp = -10°, 0°, 10°, 20°, ij-Q 0 ) on the trapezoidal wing 
measurements of the NACA (1) (A = 6; Z = 2; cp = 0°, 15°, 50°) 
on a fairly recent series of observations of trapezoidal 
wings with -55° < cp < l»-5 °, by Luetgebrune (6) (7) and on 
an unpublished DVL measurement of a trapezoidal wing 
with cp = 0° and cp = 55° • Regarding the Multhopp . 

dc fi 

thesis that -—.a, is not influenced by sweepback, this 
da 

is certainly true’ for the interval of small sweep angles 
up to 15 °, or more precisely it is not detectable since 
dc Q 

changes of — a. by 2 percent or 5 percent which would 
da 

be expected from the L-method for cp = 15 ° are hardly 
perceptible experimentally, especially- since the c a (a) 
curve Is generally not entirely straight? but this 
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assertion is also made by Luetgebrune (6) for large sweep 
angles on the basis of his measurements. In my opinion, 
however, these measurements cannot be used to decide -the 
matter, since the dispersion (explained by inaccuracies- of 
the model) of the measured c a ’ values about the mean, 
which was supposed td be independent of cp was • • 

almost ilO percent (see fig. 9 '6))> so that the ‘errors 
here were about as great as the difference in question. 
Similarly the VDT measurement (1) cannot be used, . 
particularly because the corrections to that cannot be 
checked, • If one attempts to Interpolate the Blenk 
values of o a (a) by a straight line, one gets for cp = 30° 

a decrease of c a * by about 6 percent compared to cp '•= 0°, 

while from table V one can take 10,7 percent. The DVL 
measurement (A ~ 5*7> 2 ~ 1.9) shows for cp = 35°. & 
decrease of about 8 percent, while the theoretical approx- 
imate formula would give about lij_ percent. If one were 
to draw a final conclusion from these two results. It 

do 

would have to be said: & is decreased by sweepback, 

' da. _ 

but the decrease appears to be only about 60 percent of 
the theoretical value. (When swept-back wings are. con- 
structed in the usual manner according to which the 

ci 

profile; and In particular its percentage thickness ^ 

is given in the direction of the wind, then c a ’ really 
should have another, correction because the profile sections 
in. the- direction of the. effective flow incidence - that 
is, perpendicular to. the one-quarter chord line - have a 
different thickness than- the prescribed profile.) According 
to RIngleb (12) c a - when -sweep angle Is used changes 

by the factor k = cos <p + O.723 £ (1 - cos cp ) where 

the first term is due to the changed thickness. . Since 
this influence .Is already contemplated in our calculation, 
c a ! should be multiplied . by the factor 

1 + 0,723 ^Qos^ '' <p ” * For. a thickness of 12 percent 

this means for cp - 30° an increase of c a * by 

1*3 percent and an increase of 3*6 percent for cp = h&° • 

While the above discussed question of the total lift 
is perhaps practically not so important, the practical 
engineer is especially interested in the form of the lift 
distribution on swept-back wings, especially because of its 



22 


NACA TM No. 1120 


importance for the position of the neutral point, and for 
the behavior of the wing with respect to flow separation. 
From the experimental side this question may only be 
completely answered by pressure-distribution measurements. 

To date, these have only been made by Luetgebrune who 
measured a trapezoidal wing with <p = 0° and with 
cp = ±35°. The lift distributions obtained do not reveal 
any notable influence of sweepback at all. (This fact, 
however, may be due to' the circumstance that the measure- 
ment was carried out on a wing-half with end plate at 
the center of the wing, so that the behavior at the center 
of the wing, precisely where the greatest effect of 
sweepback is to be expected theoretically, might have been 
falsified;) 

From the balance measurements only one intergrated 
value is to be taken, namely the position of- the neutral 
point on the wing chord, or for cp / 0 the spanwise 
position of the center of lift of a wing-half. In this 
way the general impression, based on experiments, that 
the neutral point may be determined from the lift 
distribution of the unswept-back wing If cp < 15 °, is 

confirmed by theory, since a deviation of l/2 percent 
of the average wing chord would be difficult to detect 
experimentally. For cp = 30 ° Kuhle ( 5 ) found In the 
NACA measurements a difference of 17 percent of the average 
wing chord between calculation and experiment; this 
enormous difference is due to a mistake in calculations, 
but the Multhopp evaluation of this measurement does show 
a comparatively great influence of. sweepback on the 
position of the center of lift, which for the most part 
agrees quite well with the results of this calculation, 
however, one must' bear in mind that for the experimental 
determination of the center of lift (in Its dependence on ), 
it is not the distance of the measured neutral point from 
the quarter-chord line (which is its theoretical position 
at cp =0) but rather its distance from the neutral point 
measured at cp a Q which is of consequence. Then one gets, 

using ^ = 3 » from table VI ( 1 ) the following values for 

the position of the center of lift: 

cp = i 5 °. a = Q« . 0£ | + . P . «?5 . g . = 0 .163 

3 tan cp “ 

cp = 30°: a = + °«775 = 0.I60 • 

3 tan cp 
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/That is, the center of lift according, to. ‘this, , actually 
is displaced some'Wftalr toward the center of the wing, as 
the.. 3 wee check increases. The. Blenk wings showed the 
same' behavior: no increase in- .a ..in going from . 

cp =.- 15 ° to cp = * 0 °, rather, a small decrease. If one 
assumes that -the experimental value of a determined for 
■<p = 15° is also approximately correct for -9 - which 
• may very well be the case,.- there are, especially for 
cp = 0 °, great differences between theory .and experiment. 




: — . — ; 

..Theory i 

Exper.imeji.tL 

a at 

cp = 0° 

L-method 

Multhopp 

Rectangle A 

= 5. : 

0.439 -. 

0 .l£l 

0.475 

Trapezoid 

= 6 , z = 2 

' ■ -k?-k 

• . 450 '. 

M 3 


As a matter ‘of fact, - pressur e-distribution. measurements , 
show that the actual lift distribution hear, the wing tips.' 
is -greater* - , than the theoretical, which may be' .explained 
by the’ influence of tip vortices Hansen (14.) gives no 
‘'neutral- point positions, pfobsbiy on account of' the. 

4 rather, nonlinear variation of c m (°a)* It mhyy however, 

be deterrfiihed here also that the Multhopp- method for' - 
cp =..lj.O° .gives too great a displacement..- - 

In summary it may.be said that the Multhopp method 
gives too large value, s for the displacement of the neutral 
point, and the sairfe-drs probably aiso true for the. 

F-method and the L-method, since experiments up to 
9 ;= 30° - apparently show no . influence of sweepback whatever 
’on’ the conter-of-lift position. Presumably this behavior 
.may be explained by assuming that. the - boundary’ layer, 
following the pressure gradients of the : swept-back wing, 
flows from the middle toward the wing tips where it is 
piled up by the oppositely directed flow about* the wing 
tips. The presence of such flow in the boundary layer 
may be clearly seen from flow pictures {tufts, and 
coloring matter in the water tunnel) made by H'an&en. 
Another indication that the differences between theory 
and experiment are to be sought in boundary- layer 
infr^uehces is the fact that the influence of- the Reynold’s 
number on .c a (a), c^Ccg), and °a max ' is . considerably 

greater- £or swept-back wings than for straight wings. 
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Finally, a. few remarks about wings without sweepback. 
With reference to the form of the lift distribution, the 
new methods may be expected to give not better agreement 
with measured lift distributions, but rather a somewhat 
poorer agreement; for while the difference between the 
old and the new distributions for large aspect ratio 
consists in a decrease of the lift, especially at the wing 
tips, the experimental distributions at the tips already 
give .values too large (compare the note on p, 22). Since 
this effect clearly has the character of nonlinearity with 
respect to a it will not be comprehended by any linear 
theory of the lifting surface, no matter how accurate. 

Hie difference between the theoretical and 
experimental values of c,* has hitherto been explained 
by saying that, even for the wing with infinite aspect, 
ratio, c a T could be expected to be less than the 

theoretical value. This hypothesis is only conditionally 
true, as is shown by a glance at figure 15. *n this 
figure the calculated values of c a » according to the 
F- method, and by the Prandtl lifting-line theory using 
c' &C B~ 2tr are plotted versus the aspect ratio A, 
together with some' experimental values obtained with very 
thin wings, which agree very well with the c» a p curve* 

According to this, there is no boundary- layer effect 
present for very thin profiles, c * a . is equal to the 

theoretical value 2tt, and the differences • obtained 
hitherto are due solely to failure to take into account 
the influence of the stirfac© effects. Naturally even here, 
especially for small aspect ratios, a perfect agreement . 
between experiment and lifting-surface theory is not to 
be expected on account of the tip effects . In addition, 
for normal wings, there is the influence of profile 
thickness, which according to the plane theory should have 
as a result an increase of c’ aoo while experiment shows 
a decrease. It Is this effect which may probably be 
correctly attributed to the boundary layer. 

In this connection a method for calculating ' c m 

similar to the F-method would appear to be desirable, 
since the relation between c m and c a can be tested 
directly by experiment. Since the induced velocity on 
the surface Increases from the leading edge toward the 
rear, a displacement of the neutral point from the one- 
quarter chord line forward is to be expected with decreasing 
aspect ratio, and this is in harmony with experience. 
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IV. NOTE .ON* THE, SWEPT- BACK WING IN YAW 


The formula • for ' the rolling 
incorrectly in almost all of the 
argument '-’being i the lift: of a s 
in yaw by ' cos2p from which the 
of a swept -back wing- In yaw is 
advancing wing half, and cos2( 
half. If A Is the total lift,- 

center of lift -of a wing half at 
due to sideslip will be 


moment due to yaw is given 
pertinent literature, the 
traight wing is decreased" 
change for the two halves 
cos2(cp - p) for the 
cp + p) for the lagging 
and If we assume, the 




the rolling moment 


L = 2 ^ £oos 2( cp - p) - .cos 2 ("cp + p)J = sin 2cp sin 2p 

. The error in this, derivation lies in the false analogy 
between the straight wring in yaw, and the swept-back wing'. 
If a straight wing is put in yaw, both the normal and 
tangential components are decreased by cos p so that 
the lift takes on a factor of cos2p. But If a straight 
wing is given sweepback, only -the effective tangential 
component Is changed by cos cp while the normal component 
remains unchanged, .so that the lift Is only to be 
multiplied by cos cp . If now the swept-back wing is 
yawed, the normal components on both wing-^halves change 
b 7. cosp k 1 and the' tangential components take on the 
factors cos( cp- p). and cos(cp+ p) • Accordingly the 
correct rolling moment .due to. sideslip is 


L = cos(cp - p) - cos ( (p. + pjj = ^ s'Incp sinp 


C T = — ^ — = — sincp sinp 

L I^pI 2 


(59) 

(iiQ) 


while the old formula gives t.wice this amount, ^fracing 
down this error in the literature Is further complicated 
by the circumstance that the rolling moment is not always made 

dimensionless with Si as in this report, but with 
instead. ^ 


Fb 
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So for instance, in the much used summary by Schlichting ( 13 ) 
there is a formula for for small cp and (J which 

is formally the same result as in equation (ij.0). Never- 
theless the Schlichting formula contains the above error 
since was used in nondimens ionalizing. 

The only work known to me which avoid the above 
fallacy is Multhopp’ s paper (10), but since the error, 
even here, is not expressly emphasized, and rolling moment 
due to sideslip is only given in the form of an integral 
to be evaluated in every individual case over a lift 
distribution which has to be calculated for this case, the 
Multhopp results have been less noticed than some others. 

In addition, Multhopp is concerned with the oblique 
position of the free vortex sheet, which causes an 
additional angle of attack distribution Aa - a^(3 tancp 

where cu is the induced angle of attack distribution 
for nonoblique flow incidence j the positive sign is to 
be taken for the advancing wing, and the negative for the 
retarded wing. While now, according to Multhopp, the 
lift distribution belonging to Aa must be calculated, 
and from' it the corresponding rolling moment due to side- 
slip, one may by a somewhat cruder procedure obtain also 
in this case a closed formula. 


To this end we assume the lift distribution of the 
wing without yaw to be approximately elliptical. Then 
a ^ is constant along the span, and equal to 

C„ 

likewise the absolute value of Aa is constant along 

the whole span; but the sign changes at the center from 
plus to minus. The lift distribution corresponding to a 
discontinuous distribution of a like this, is, for a 
wing half, approximately like the usual lift distribution 
of a wing of half the span. That Is, the total lift 
coefficient associated with the Aa of a wing half, if 
we use the elliptical conversion formula, and c ’ aco = 2ir 
is given by ** 


2tt coscp 


A + 
2 + 


Aa = ±rr(3 sincp af= ip sincp 


2C ( 


|+2 


A + ij. 


(hi) 



naoa m mb, iiao 


2? 


From this we get for the coefficient of the corresponding 
rolling moment: 

i3sln ' p rf-£' 

Together with the moment (ij.0) of the bound vortices, this 
gives for the total rolling moment due to yaw resulting 
from sweepback the convenient formula 

5 T =C *( 0 * 5 + rfaW • ^ 


V. SUMMARY 


Two methods have been developed for calculating . 
the lift distribution over the span, which takes into 
account the”' influence of the distribution of the 
circulation- over an area better than the Prandtl lifting- 
line theory, and which may both be used on wings with 
sweepback* For swept-back *vings, the first method, 
called the F-meth’od, is numerically very laborious, and 
therefore, serves only as a check of the simpler 
L-method. The check is very good, even at large sweep 
angles* The computations for a straight wing require 
about 3 hours, for the swept-back wing, about 8 hours* 

A series of examples .was calculated numerically, from 
which the following conclusions may be drawn; 

1* For the straight wing the new method gives a 
noticeable decrease of c* when compared with the 

ci 

Prandtl method, which for example at A = 5 amounts 
to about 8 percent. If the profile is very thin, the 
experimental values seem to lie very well on the new 
c* a curve. According to this, the difference between 

the old lifting-line theory and experiment, in the case 
of very thin profiles, is not to be attributed to a 
decrease of the theoretical c» s 2if because of 

- «-Co 

boundary layer, but comes from neglecting the two<* 
dimensional distribution of the circulation,*- 
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2. In figure if the- ratio- 


°V? 

° ’a P 


of the new c 1 . 


to the old is plotted versus the aspect ratio A . Since 
the curve is pretty flat in the interval of normal aspect 
ratios, the conversion formulas for changing from one 
aspect ratio to another are not essentially changed. She 
curve may, however, be used for correcting the rolling 
moments caused by antisymmetric al angle of attack 
distributions. For this purpose, one would have to 
multiply the moment calculated by the old theory by the 


C5 * K 

value of — - — soaled off of- the curve at the point *3 l 

c’a p ' 2 

Thus for example one would get for the rolling moment due 
to sideslip resulting from dihedral a decrease of 
15 percent while measurement gives 18 percent. This Is 
for a rectangular wing with' aspect ratio of 5* 


5 , For swept-back wings a comparison of the L -method 
with the method' of Multhopp (10) was carried out. The 
L-method gives a greater decrease of c ' a and a smaller 

displacement of the neutral point caused by sweepback 
than the Multhopp method; moreover. It converges more 
rapidly. 

4* Available experiments show less change of c f 

caused by sweepback than predicted by the L-method, and 
either no displacement at all or very small displacement 
of the center of load in. the spanwise direction. These 
differences between theory and experiment are presumably 
to be explained by a movement of boundary-layer material 
toward the wing tips, and this is verified by flow 
observations . 


5 . The formula customarily given for the rolling 
moment due to sideslip of the bound vortices of a swept- 
back wing rests on a fallacy, and gives values which are 
100 percent too larg 6 . In section IV a convenient 
closed formula is derived for the rolling moment due to 
yaw resulting from sweepback, iri which formula the 
influence of the free vortices is also taken into account 


Translation by H. R. Grummann 
McDonnell Aircraft Corporation 
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SUPPLEMENT 


In comparing my evaluation of the VET measurement 
with the Multhopp evaluation ((10) fig. 3)> 1 could not 
find the point corresponding to the wing 24.-15-0. 

Writing to Mr, Multhopp about this, I was informed that 
he had omitted this point because he did not think it was 
correct. At the same time he sent me the accompanying 
amplified diagram in which there Is, In addition, a new 
point 1, 

This figure seems at first to verify the correctness 
of the Multhopp calculation, this al30 if one only 
considers the inclination of the theoretical straight 
line. In my report only wings 5 and 7 were used, because 
these were the only ones for which nothing was changed 
but the angle of sweepback in starting from cp =0°; 
for all ether wings the profile and the twist were changed 
at the same time, for no. 1 even the plan form was 
changed. Ihe thesis that sweepback causes no essential 
displacement of center of lift is thus confirmed by the 
only comparable measurements 5 and 7 (and by Blenk (3))« 
Against this conclusion may be adduced the other measure- 
ments, as well as the fact that by extrapolation to 
9=0° one obtains centers ox lift located very far out. 
(The allusion in my report to tip effect may not be 
sufficient as an explanation as Multhopp correctly 
remarks.) Finally, if one considers that the calculation 
of a from measurements of c m rests on the assumption 

that the neutral points of the individual profile sections 
are not changed by sweepback, one must probably say in 
conclusion that a final answer to the question cannot 
yet be given on the basis of the experimental results 
available at the present time. 
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TABLE I 





n 

■■ 

1 

n = 3 

n = k 

0 

0 

■ 0 . 

. 2.6151 

*1.1*1 i|2 

1.062k 

-0.5000 

■M 

- 

1 

4.5S89 . 

-2.1053 s 

1.5663 

-.7193 

9m 

1 

2 

•2 ; 3Skk 

-.2363 •’ 

+.0751 

-.0253 

«T 


z 

•6575 

. 2.001*6 

-i.5k02 

.7109 

1 

• 2 


-1.1*11*2 
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CONTINUATION' OF TABLE III 
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TABLE V 
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The coefficient c» as a function of the swept-back 

8 . 

angle (according to the L-method) . • 



TABLE VI 


IMIMI 

Mlth MnlthooM b qua rim fonnula 

Graphic *1 

Method 



F-method 

__I 

L-W 

»thod 


L-b 

nethod 


Wing 

^0 

<? 

■1 

Aft 

£ Aft tan y 

a 

Aft 

m 

D 

AA 

£ &a tan f 

B 

Aft 

£ H* tany 


0 

0.451 

0 

0 

0.440 

0 

0 

0.439 

0 

0 

0.438 

0 

0 

A = 5 

15 




.451 

.on 

• 007 

.450 

•on 

.007 





50 







•463 

.024 

.035 




Z = 1 

45 

.505 

.054 

.155 

.403 

.043 

.107 

.461 

.042 

•105 

.472 

.054 

.085 


0 

.429 

0 

0 . 




.424 

0 

0 

.423 

0 

0 

A- 5 

15 







*433 

•009 

.006 





50 







.443 

.019 

.027 

.440 

.017 

.025 

2 “ 2 

45 

.4fl5 

.056 

.140 




•457 

.033 

.082 

.450 

.027 

.068 

A = 6 

0 

.450 






.424 






2 = 2 

50 







.447 

•023 

.O40 




Ac 10 

0 

.451 






.427 






2 = 2 

45 







.472 

.045 

.225 





NATIONAL ADVISORY 
COMMITTEE R* AERONAUTICS 


Bio distance a from the plane of symmetry of the lift center 
of gravity of ft wing half In terms of aemlftpan#- 


co 

oo 


NACA TM No. 1120 












Figure 1. The vortex system of the lifting surface 

(a) Straight-rectangular wing 

(b) Swept-back rectangular wing 
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Fig. 3 



•Figure 3. Lift distribution of the 3 rectangular wings = 1, 5, 10, 
Comparison of Prandtl's theory of the supporting line with the 
lifting surface method. 




^ap ^ap 

ratios 0-^“10 ( C a ^ = lift coefficient according' to the F-method, 


C a ^ - lift coefficient according to Prandtl's theory of the 
supporting line with C_ = 29f)« 


Fig. 4 NACA TM No. 1120 
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Figs. 5,6 



Figure 5. The 'lift distribution G (Y) of the 45° rectangular 
swept-back wing-A = 5 as the sum of. 2 distriburions (} (Y) 
and C (Y)V 1- Y z (F-method). 



Figure 6. Lift distributions of the rectangular wing-A. = 5 with 
the sweep-back angles s? = 0°, 15° and 45° according to the 
lifting surface and supporting line method. 
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Figure 7. The vortex model of the L-method 

(a) Wing withbut sweep back 

(b) Swept-back wing 
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Figs. 9,10 
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Figure 9. Lift distributions of the rectangular win g«/\_ = 5 with 
the swept-back angles v = 0°, 15°, 30° and 45°, according to 
the L-method. Dashed line; result of Multhopp's method for 
(K = 1, m = 7 ) for = 0° and 45°. 



Figure 10. Lift distributions of the trapezoidal wings, -A- = 5> 

Z = 2 with the swept-back angles - 0°, 15°, 30° and 45°, 
according to the L-method. Dashed line; result of Multhopp's 
method for ; (K = 1, m = 7) for ^ = 0° and 45°. 




Fig. 13 
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theory of the supporting one 0^^ = 2"K) compared with measure- 
ments on very thin profiles. 


